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ABSTRACT 

Spacecraft measurements of core ions on L=4-7 field-lines typically show trapped ion 
distributions near the magnetic equator, and frequently indicate field-aligned ion streams at 
higher latitudes, (e.g., Olsen et al. [1987]). The nature of the transition between them may 
indicate both the microphysics of hot-cold plasma interactions and overall consequences 
for core plasma evolution. We have undertaken a statistical analysis and characterization 
of this interface and its relation to the equatorial region of the inner magnetosphere. In 
this analysis, we have characterized such features as the equatorial ion flux anisotropy, the 
penetration of field-aligned ionospheric streams into the equatorial region, the scale of the 
transition into trapped ion populations, and the transition latitude. We found that most 
transition latitudes occur within 13° of the equator. The typical values of equatorial ion 
anisotropies are consistent with bi-Maxwellian temperature ratios of T±/T\\ in the range 
of 3-5. The latitudinal scales for the edges of the trapped ion populations display a rather 
strong peak in the 2 — 3° range. We also found that there is a trend for the penetration 
ratio, the anisotropy half width, and the transition scale length to decrease with a higher 
equatorial ion anisotropy. We may interpret these features in terms of Liouville mapping of 
equatorially trapped ions and the reflection of the incoming ionospheric ion streams from 

the equatorial potential peaks associated with such trapped ions. 
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Spacecraft measurements at tprr€sshow a latitudinal transition of core ion (0-50 eV) 
fluxes in which trapped ionsj^re seen centered about the magnetic equator and counter- 
streaming "ions are seen at off-equatorial latitudes (e.g., Olsen et al. [1987]). This indicates 

thft j mor a rr mf tp a p p mg the field-aligned flowing ions in the equatorial regions, which is an 

^ •*r" 

important aspect of plasmasphere refilling* The nature of this interface is quite interesting 

4 C* *t tCi**y . 

from the standpoint of both the microphysics of hot-cold plasma interactions and overall 
consequences for core plasma evolution. 

There are a variety of core plasma observations near the equator in the vicinity of the 
plasmapause. These observations often reveal plasma pitch angle distribution with peaks 
at 90° (pancake or trapped distributions) and characteristic energies of several electron 
volt$(e •g. Horwitz and Chappell, 1979; Comfort and Horwitz, 1981; Horwitz et al, 1981b; 
Olsen, 1981; Olsen et al.,1987, Sagawa et al. 1987). Since these ions originate in the 
ionosphere, it is clear that they have been subjected to perpendicular heating wMchwe, p 
presumably responsible for the waves associated with coexisting ring current populations. & 


^Sagawa et al. [1987] find the occurrence frequencies for such trapped distributions observed 
[^within 5° of the magnetic equator to be in the T=3.5-5 range. They also find broadened 
field-aligned distributions and conics middle latitudejin the outer plasmasphere, 

which may indicate wave-particle interactions at lower altitudes along these field lines. 

Singh and Hwang [1987] have examined the equatorial trapped ions in the refilling 
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processes and suggested that the equatorial heating is likely to set upidownward electric 

' 

field due to the difference pitch angle anisotropies of the trapped ions and electrons at 
the equator. The electric fields and the associated potentia\larf"lhe equator were estimated 
to be large enough to reflect the field-aligned flows reaching the equator (Singh and Torr, 
1988). Lin et al. [1994] have used Generalized SemiKinetic Modeling (GSK) to simulate 
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development of the potential structure and the reflection of the incoming ionospheric flows, 
and demonstrated effects of hemispheric decoupling and asymmetric refilling owing to such 
equatorial potential peaks. 


In recenjt years, several new characteristics of thermal plasma distributions have been 

US'*/ 1 ) 

observed Jjy'DE-l /RIMS (Dynamics Explorer-1/ Retarding Ion Mass Spectrometer) [e.g. 

Chappell et al,1981]. The RIMS instrument je^designed to measure the energy and angular 
distributions of low-energy ions for various ion species. In particular, analysis of RIMS ion 
data provides estimates of wide ranges of thermal ion bulk parameters, including the ion 
density ( 0.1-10 6 ions/cm 3 ), temperature (0-45 eV), and bulk flow (> 0.5 km/sec) in the 
inner plasmasphere and ionosphere and the^sp^rffic^ktn pitch angle and energy spectral 
characteristics in the outer plasmasphere and plasma trough for a mass range 1-32 amu. 

The RIMS instrument consists of a central electronics assembly And three separate sensor 

heads, one mounted viewing perpendicular to the*spin axis^and two each mounted with . . 

' * 

field of view parallel and anti-parallel to the gj^aeeCfaft spin axis. ^ 

^ k&ydl ) 

An unique opportunity offered by DE 1 core ion observations was that during certain 
periods the satellite orbit essentially skimmed a large latitudinal segment of approximately 
constant L shell of ~4.6 (in the vicinity of plasmapause) about the equator. This allowed 
excellent measurements of the latitudinal profile of the particles and fields in this region. 

Olsen et al. [1987] utilized such orbits to examine the characteristics of the low energy ions 
and waves in this region. They observ^trongly confined light ions about the equator and 
equatorially confined broadband, low frequency electrostatic waves. Olsen et al. [1987] 
also show cfata in which there is either a gradual fade in the intensity of field-aligned flows 
or what appears to be field-aligned ion beams reflecting from a thin potential structure. 

The low penetration ratios of field- aligned streams into the equatorial regions in some 


cases imp 



ctrostatic potential layers which are indicative of important laminar 


hot-cold plasma interactions. Such layers also should lead to hemispheric decoupling and 
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hence have major consequences for the plasma evolution along these refilling flux tubes. 
As noted above, Lin et al. [1994] have simulated core plasma evolution with equatorial ion 
and electron heating. They found that, around the equator, electrostatic potential barriers 
developed and resulted in reflection of incoming ionospheric flows and hence strongly influ- 
enced the bulk parameter profile evolution. Olsen et al. [1987] have obtained the electric 
field of about 0.1 fiW / m by referring to the model of Liouville mapping and bi-Maxwellian 
distribution for equatorial trapped ions (Olsen [1981]). Such potentials have been proposed 
as an explanation for the changes in field- aligned ion behavior in the DE 1 data set. 

In this report, we present the results of a statistical study of the properties of the 
transition between the trapped and field-aligned ion fluxes in the equatorial region, using 
DE 1/RIMS data of mid-1982 and late 1983, orbit was particularly suited to lati- 

tudinal profiles with nearly constant £, with apogee near the equator. In this analysis, we 
have characterized such features as the equatorial ion flux anisotropy, the penetration of 
field-aligned ionospheric streams into the equatorial region, the scale of the transition into 
trapped ion populations, and the transition latitude. Statistical properties and correlations 
of these parameters will be described. These features will be interpreted in terms of Liou- 
ville mapping of equatorially trapped ions and the reflection of the incoming ionospheric 
ion streams from the equatorial potential peaks. 

In section 2, we will describe data analysis methods used including criteria for selected 
DE-1 measurements and parameter definitions. In this section, the method of Liouville 
mapping will also be briefly described. Section 3 describes the results including statistical 
occurrence trends and corrections between defined parameters. A discussion and conclusion 
are given in section 4. 
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2. METHOD 


2.1 Criteria for selected DE-1 measurements 


V'oy 
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A statistical study of the properties of the latitudinal transition between the trapped 
and field-aligned ion fluxes in the equatorial region has been investigated by utilizing RIMS 
core ion data from DE-1 orbits which nearly ”skim” along the L=4.6 shell. The criteria 
for selection of appropriate passes are first that they involve a transition from^trajrj^ed 
counter-streaming H + ions and that the variation of L-shell near this transition is small. 

Figurejl and 2 show two typical spectrograms for DE-1 passes through the equatorial 
region when the spacecraft apogee was near the equator. These two panels are spin-time 
spectrograms of ion count rates data taken at 0 V retarding potential. The center line is 

Ate* 

the ram direction which is nearl y^ero degre es p itch angle^ The minimum pitch angle (near 
0° ) is indicted by white line running along the center of the spectrogram. The white line 
near the bottom of the panel indicates the spin phase for maximum pitch angle (near 180°). 
The i^^es ^count rates^are color coded with dark blue representing low count rates and 
red representing high count rate. The color bar on the left hand side shows this scale. Note 
that the two spectrograms have been scaled differently to emphasize the desired features. 
The equatorially trapped proton distribution appears as the red region between the two 
white lines in both figures, indicating a flux peak^^O 0 pitch angle. Prior to the spacecraft 
encounter with these trapped ions, the proton distribution is primarily field-aligned (fluxes 
peaked along the white lines) in both cases. The principal difference between the two cases 
is that Figure 1 shows an abrupt transition from field-aligned to trapped protons and very 
low penetration of the field-aligned protons into the trapped ion region while Fig 2 shows 
a more diffuse transition and higher penetration of field-aligned protons into the trapped 
ion region. It is evident that analysis of such transitions in this region will aid significantly 
in understanding the nature of the large-scale core plasma evolution on these flux tubes. 



Regarding our criterion for the acceptability of DE-1 passes for this study, the unique 
capability for sampling the distributions of core ions along the field line in the vicinity of the 
magnetic equator is illustrated for two example orbits in Figures 3 and 4. Figure 3 shows 
the DE-1 orbit for the period 09:00-12:00 of May 6, 1982 (82- 
superposed on model dipole field lines projected to a geocentric distance-magnetic latitude 
plane. The orbits, of course, do not stay precisely on an L=4.6 shell, and this aspect 
does introduce the possibility of L-shell variations, as opposed to the desired latitudinal 
variations alone, affecting the observed plasma distributions. Therefore, we calculate the 
L-shell variation between A = +10° — 0° or A = —10° — 0° (depending on which hemisphere 
the transition point is located), which generally brackets both the equator and the observed 
field-aligned/trapped transition locations. 



LL — |Lx=±io<> — La=o°| (1) 

For the cases in Figures 3 and 4, the LL values for 82-126 and 83-295 were 0.13 and 

HO 

0.36, respectively. ^ cases with LL values greater than 0.36 were considered. 

The observation periods for this study were May to September, 1982 and September 
to December, 1983. Over these two periods, 43 passes satisfied our LL criterion for nearly- 
skimming orbits. 

2.2 Extracted parameters for survey 

In this survey, we have obtained five parameters for each pass for characterization of 

this region. These parameters are the equatorial ion anisotropy, anisotropy half-width, the 

transition latitude, the transition latitudinal scale length, and the penetration ratio for 

incoming field-aligned streams entering the trapped ion region. Figure 5 shows plots of 

count rates at different spin/pitch angles and their ratios together with parameter defini- 

tions for the 82-126 data. The data we have fcensaysd'is actually organized in terms of spin 
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angle, while we need the count rates in terms of pitch angle, and the relationship between 
spin angle and pitch angle changes somewhat during the equatorial pass. Therefore we 
have used the following approximate approach. 

For the maximum anisotropy and latitudinal width parameters, we used and plotted 
(top panel of Figure 5) the ratio of the count rates at spin angles corresponding to pitch 
angles of 90° and 45° for the spin- pitch angle relationship at the magnetic equator indicated 
from the orbital data. For example, if the white line which represents the minimum pitch 
angle is located at spin angle 15° while it passes the magnetic equator, the trapped ion 
population is taken in the range of spin angle 6 = 105 ± 20°, where 105° is the 90° pitch 
angle for trapped ions in addition to spin- pitch angle difference 15°. 

To estimate the field- aligned count rate variations along the orbit, we have used and 
plotted (middle panel of Figure 5) the count rates which corresponc^o the minimum pitch 
angle at the transition location. This transition location is initially determined roughly 
from the microfiche spectrograms as the apparent edge of the trapped ion distribution 
region. 

To obtain the actual trapped ion transition location, we have used and plotted (bottom 
panel of Figure 5) the ratio of the count rates at spin angles corresponding to pitch angles of 
90° and 45° for the spin-pitch angle relationship at the spectrogram determined transition 
location. 

For cases where the orbit follows the L-shell closely, as in Figure 5, the spin-pitch angle 

relationship is nearly constant, and thus the difference in the ratio in the top and bottom 

panel is slight; but in other cases it can be significant and it is necessary to do these 

two different calculations (both the transition point and the equator). For the parameters 

determined from the data of the type in the top and middle plots in Figure 5, we have 

performed a three-point weighted average smoothing (solid lines) of th^fninute samples 

(dashed lines) before computing the parameters as below. 
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These parameters are defined as follows: 


jy 


Equatorial Ion Anisotropy 


/ 


According to Lin et al. [1994] and earlier modeling efforts, it is known that the equato- 
rial ion temperature anisotropy will affect refilling and interhemispheric plasma exchange. 
Hence, the equatorial anisotropy will be a determining factor in this study. Here we define 
the maximum equatorial anisotropy A max from the variation of countrates as: 


A-..; 


A 

s*-ma.x — 


CR(a = 90°) 


( 2 ) 


CR(a = 45°) 

where CR is the ion count rate and the ratio is the maximum anisotropy along the orbit. 
The location of A max is generally near or at the nominal equator. Note that CR(a = 45°) 
instead of CR(a = 0°) is chosen to be the denominator here. This is because, for example, 
when both counterstreaming and trapped populations have similar count rates, the value 
of be do se t° 1? while the distribution may not be isotropic. 

We note here that A max is not always located at the nominal equator obtained from DE 




orbityalfitude data. In this study we have tried another method to determine the location 
of tne equator. Considering a dipole field, we mapped the measured B magnitude to a 
consistent L-shell value at the corresponding latitude. In other words, we corrected the 
measured B field to compensate the variation in spacecraft radial distance along its orbit 
to better determine the time of equatorial crossing based on corrected B field strength 
minimum. We found that neither the location of B m i n nor the location of the magnetic 
equator obtained from the simple model used in the DE orbit data is always consistent 
with the location of A max . However, we normally regard A max to be the physically most 
meaningful parameter (especially in relation to the field- aligned penetration ratio), and so 

we have chosen its location to be the effective magnetic equator in this study. 
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Anisotropy half width ( W ) 


To estimate the rapidity of the temperature anisotropy decline from the maximum 
anisotropy location, we define the Anisotropy Half Width (W) as the latitudinal distance 
between the locations of the maximum anisotropy value and the half maximum anisotropy 
value. 

Transition latitude for trapped ions ( X tT ) 

The Transition Latitude (A tr ) is defined as the latitude of the first data point toward 
the equator of three consecutive data points which have anisotropy greater than 1.2. This 
location is intended to indicate the entry into the edge of the trapped ion region. The 
value A > 1.2 is chosen by making a comparison between the spectrogram from microfiche 
and the ion anisotropy profile. The requirement for three consecutive points with A > 1.2 
is to discard potentially noisy or intermittent excursions above the 1.2 threshold. 

The transition latitudes obtained for the passes of 82-126 and 84-068 (in Figures 1 and 
2) were at -6.07° and 12.5°, respectively. 

Transition Latitudinal Scale Length (X,(trapped)) 

To indicate the sharpness of the trapped ion transition, we define A ,(trapped), the 
latitudinal scale length associated with the trapped flux decrease at the transition, as 

^ ^ m * 

dCR ft ^ 

• X, (trapped) = CR tr an, x ^ **’ ( 3 ) 

where is evaluated about the transition latitude by taking the slope of CR from 
A tr — 1 to A tr + 1 and C R tra na is defined by the average count rate in the same range. 

The latitudinal scale lengths of the above passes on 82-126 and 84-068 are 1.29° and 


5.63°, respectively. 
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Penetration Ratio for Incoming Field- Aligned Ioris (P) 


A measure of how much equatorial penetration of the field-aligned ions occurs reveals 
information on the electrostatic potential structure along the field fine. We define P, 
the penetration ratio of counterstreaming ions, as the ratio of the minimum field-aligned 
flux observed between the transition latitude and the equator to the field- aligned flux ^ 
observed just outside the transition location. To define P, we decided that the minimum 
flux observed between the transition latitude and the equator would be a better choice for_ 
the numerator than the field-aligned flux precisely at the equator or the A max location. , ^ 
This is in part because the largest potential barriers which will reflect the incoming flux 
could in principle be at slightly off-equatorial latitudes. Hence, the minimum flux between 
the transition latitude and the equator could be correlated to this largest potential barrier. 


Obtained parameters for the example passes 


A 


\ 




The penetration ratios for the passes of 82-126 and 84-068 we ^fh2g ^tnd 0.78, respectively. 

V£.< es 

SO 

Table 1 fists the obtained parameters for the two example passes discussed here. 



2.3 Assumptions for the relation between A max and the equatorial temperature anisotropy 
As mentioned before, the equatorial ion anisotropy A max is a determining factor in this 
study. Here, we will introduce an approximate method to find the relationship between 
A max as defined here and Tx/T\\(eq) for a <da toote 4» b i-Maxwelfian ion distribution. For an 


aperture 



arbitrary ion distribution function f{y ), w e can calcul a t e the count rate into the detector 


l, 




CR(a) = Co| J df> J sin9d.0v 2 dvf(v)v • d\ 

11 A 


( 4 ) 


10 



\r4r** * ■ * ^ d tf-c? art^oo/’ r«d'p r »r'Vir > d-- 

I (T» 4 **' i A*» % /^ ^ v^/ ^#»v ^ 

L €-eJL *, /et> £€*4. 


d G'b art^Oos f‘^d\*»r - Vir > £-' (*s>'f 4 

J > « * .. 


Table 1: Parameters of two DE-1 example passes 


Parameter 

82-126 

84-068 

V 

1.29 

5.63 

p 

0.28 

0.78 

A tr 

-6.07 

12.5 

A 

•f*- max 

10.12 

2.99 


^,k Ai 4*> tL^t***^*' dtt^c. 

where 0 , <£ and d are^ie polar angle relative to the magnetic field, the azimuthal 
angle and the unit vector n er mol to the LtetElLm 1 , respectively, «a nd -u is the ion’s velocity 

vecto^The angular integration is over the range subtended by the detector aperture fan. 

^ ft 6/ AQt &'»'<— 

(7 0 is a constant determined by i ho - r e lat ie nchi^ -b otween iun fluj p m f d inn nwwnt pnts af The 

polar angle 0 is the angle between the magnetic field direction and the normal vector of 

the detector. ^ 

ri V 

For RIMS, the ^AglUacj^pe«pi«Nifie solid angle decreases with increasing ion energy £'**' 1 
(Chappell et al. [1981]). .W f he b e am c n ugj urnF mass uurnhii 1 'iiieie gse , t he angular 1 


Boardsen (private communication, 1993) has suggested the following approximation for 
the energy dependence of the RIMS % txteJytafciVa solid angle: 


S = So- exp(-E/kT d ) . * 

Where E is the energy of the ion and kT d it> lliL Lfllllli^ uiug'^TrHwe-drfcerfcw (about 60 
eV for RIMS). U * + +<• C Kfoiv"*. 
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For present purposes, we will assume the effective MAofwF solid angle is a function of 
the ion’s energy ( S = S(E)) as indicated above. For qualitative purposes here, we will 
further assume that the detector solid angle is small enough that the integration can 
be approximated as a simple constant multiplication, giving 


1 

l 

j\* 

fJ 

L **: 


- * 


* Vi 

<5 




CR(a) = Co -So- 1 v 3 dvf(v)exp(—(E - q$ p )/kT d ) 

where $ p is the spacecraft potential, q is the ion charge, and a is the pitch angle. 
We will further assume that the ion distribution is a bi-Maxwellian distribution: 


( 6 ) 


x (m/2A:7r) 3 / 2 
f ( v ) = n ° exp{-{ 

v 2 ii T j - 


cos 2 a sin 2 a 


+ 


)E] 


( 7 ) 


After substituting the ion distribution function into equation (6), we can integrate with 
respect to energy from E m in to oo and obtain the count rate as: 


CR(a) = C 1 exp[Emin/T d ) ■ [°° [E • exp(-GE /T ± )}dE 

"Emin 


= -Cl • Tl . exp[Emin/T d } • exp[-GE min /T ± \ • ( 


E 

J - / m\r 

~tT 


■MG 


( 8 ) 

( 9 ) 


where 


e cr a fc - gfla ; 


G=^+(^-l ) Wa + 1 

T d Tu 


(10) 


Emin H t i hf* gy 'Aj , Ci absorbs the constants 

of the distribution function and the detector. The lower limit of the energy integration 
E m i n indicates that only ions with greater energy than q$ will be detected. 

We did not try to measure the perpendicular temperature for the present analyses, 
although this is in principle possible during these periods using the Z-heads of RIMS (e.g. 


Olsen et al. [1994]). However, if we assume that E m i, 
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the spacecraft 


^gatively 
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charged or unb ias ed^then the ratio of the ion count rates at pitch angles of 90° and 45° is 
related to the ion temperature anisotropy as 


CR( 90) f£ + Krj~ + 1 ) 2 

CR( 45) 1 ^ + 1 J 


( 11 ) 


For cases where the perpendicular temperature is much less than 60 eV, ^ < 1, the 
relationship is 


CR(9 0) (rj' + l) 2 

CR( 45) “ 4 


( 12 ) 


For cases where the perpendicular temperature is near 60 eV, i. e., w 1, the 
relationship is 


C£(90) (t^ + 3 ) 2 

CR{ 45) ” 16 


(13) 


2-4 A simple model for the interpretation of equatorial core ion data 

We will also discuss in greater detail the correlations among the parameters we defined 
in previous section ( A max , A tr , A a , W, P) and possible physics that these trends indicate. 
We will do this by the aid of equations obt ained by Liouville mapping of the ion distribution 
function along the field line. 

Here we will assume a bi-Maxwellian ion distribution function at the equator, and ob- 
tain the distribution function along the the same field line. Let the equatorial distribution 
be 


r/tr Tr ^ (m/2fc7r) 3/2 

f(V\\,V x ) eq = n„ • ^ exp(- 

V r ll' 


mV I, 2 


11*9 


2kT\\, 


mV L q 

2kT ± 


) 


(14) 


II eqj-leq “• — - I|e« 

To map this distribution to an off-equatorial location A, we assume constancy of the ion 

magnetic moment and energy conservation and rearrange the resulting mapped distribution 
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function in the form of a bi- Maxwellian outside of a “forbidden region” in velocity space 
associated with the potential (e.g., Chiu and Schulz [1978]; Whipple [1977]). We can obtain 
the potential distribution along the field line, if we specify a bi-Maxwellian or isotropic 


equatorial electron distribution, require quasi-neutrality, and assume that the forbidden 
region smoothly filled in by the ions confined in the region between the eqatorial 


potential and the mirror points (e.g., Olsen et al. [1994]). 


Under these conditions, the parallel temperature is constant along the field line, while 
the variation of the perpendicular temperature is given by "ZT ~ * “ ^ 

* n# 7 r„ 



■leg 


[1 - Rb\Rt,„ + Rb 


(15) 


where Rt, 


(/» 


> 1 


The variation of the ion density along the field line can also be obtained this mapping 
(e.g. Olsen et al. [1993]). 


a. Estimation of from A r 


If we further assume that the electron temperature is isotropic and satisfies Boltzmann 


relation for the electrons, we can obtain the potential as it varies from the magnetic equalor 

^ p lc0 V 

i 






A 4* = 


k T e \\T\\i 


■ln[ 



e[T e || -I- T\\i] 1 Rb + (1 - Rb)Rt< 




(16) 



= B eq /B a and Rt = T±i/T\[i_ 


Assuming the potential distribution is dictated largely by differential anisotropy be- 
tween the ions and electrons, we can use equation (16) to relate the maximum anisotropy 


( A max ) and potential drop at a given location (A$ ). Therefore, at a given location and 
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electron temperature, we can find the relationship between A$ and Rx eq and hence A$ 
and A maw , 

b. Relationship between a reflecting potential drop (A$) and P 


We can estimate the potential drop associated with the reflection of ion streams for 
a given penetration ratio as follows. Let us assume for simplicity that the ion stream- 
ing distribution entering the interface and reflecting region is the forward portion of an 

Op*K>ArS y . 

undisplaced bi-Maxwellian distribution: r 


/(Vj,,Vl) = no 


(m/2kir) 3 / 2 

•JWl 


exp(- 


_mVl _ rnVl 
2kT„ 2kT, ’’ 11 p 


(17) 


/(V [| , Vj_ ) = 0, Vj| downward 

If the ions enter the interface with parallel velocity Vf| 0 , their energy is related to that 
inside the potential region 


V & = V l\" + 


2eA$ 


m 


where A$ is the increase in potential over that outside the interface. 
Thus, the distribution function within the potential region is: 


(18) 




2kT\\ 

From the above equation, it can be seen that only the particles with VJjq > 
reach the equator. 


(19) 


2eA^ 


can 
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The fluxes at the ionosphere and the equator can be obtained by integrating the distri- 
bution functions, giving a ratio whose maximum value we may identify with our penetration 
ratio P: 


eA$ 

Note that 7]| here is the parallel ion temperature of the the 


P = 


Flux(eq) 

Flux(incoming) 


= exp[ 


'] 

incoming flow. 


( 20 ) 
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3. RESULTS 

In this section, we shall present trends observed in the five parameters defined in 
previous section. These five parameters are: (1) the maximum anisotropy (^4 mai ); (2) the 
transition latitude (A tr ); (3) the transition latitudinal scale (A,); (4) the anisotropy half- 
width (W); and (5) the penetration ratio (P) for the incoming field-aligned ions. The trends 
presented here are intended to provide a statistical picture of the core ion distribution in 
the near-equatorial inner magnetosphere. Discussion of the results in terms of simple 
models involving potential barriers and Liouville mapping of ion distributions along the 
field line is also contained in this section. 

1. Occurrence trends of the ion parameters 

In this section, we present the occurrence histograms for the various individual pa- 
rameters. In all cases these are presented as bar graphs indicating the total number of 
occurrences within incremental bins spanning the range of the observed values of the par- 
ticular parameter. We presen t#^! he occurrence trends of the ion parameters for 43 

cases. 

Occurrences of Maximum Anisotropies 

The dominant influence on the latitudinal ion distribution 
maximum anisotropy (nominally the equatorial ion anisotropy) 
should determine, in part through Liouville mapping, the tran 
should also control the field-aligned flux penetration ratio through its influence on the 
equatorial electrostatic potential distribution. Here, we show the occurrence histogram for 
A max in Figure 6 (first panel). 

FigureG (top panel) indicates that the equatorial anisotropy ratio peaks broadly about 
->4ma* = 4, while 3 cases with A max > 10 are observed. The average value is indicated to 
be < A max >= 5.19. 
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According to equation (11), for the case that T± is much smaller than th e effective 
energy of the detector kTd (& 60eF),. we can obtain the temperature anisotropy from 


Ix/Tj| = 2 \J A max — 1 


( 21 ) 


Hence, the peak occurrence of A max = 4 would correspond to a typical equatorial T±/T\ \ 
of about 3. In the case of the RIMS instrument, the effects of the energy-dependent solid 
angle imply higher Tj./Tj| than indicated by (21). For the case that T±/Td ss 1, the average 
equatorial temperature anisotropy would be about 5. 




Penetration Ratios 


One of the most interesting parameters obtained for this region is the penetration ratio 
P, which may be important for indicating the presence of electrostatic potential barriers 
and hemispheric decoupling. Figure 6 (second panel) indicates the occurrence distribution 
of the penetration ratios P, with P ranging from 0 to 2. 

Most of the penetration ratios are less than 1, which implies the existence of positive 
potential layers causing reflection of the incoming ionospheric streams, and diminution 
of the stream inside the transition point. However, note that 10 out of 43 cases have 
penetration ratios greater than 1. This could indicate an accumulation of large field- 
aligned fluxes within an electrostatic potential well inside the transition location; such 
fluxes might not be instantaneously related to the incoming flux outside the transition 
region. The average penetration ratio for those P > 1 is 1.26 and the average value for 
P < 1 is 0.57. 


Transition Latitudinal Scales 
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The transition scale A ,(trapped) is a measure of the sharpness of the trapped ion 
latitudinal distribution near the transition location. Figure 6 (third panel) shows the 
occurrences of the transition latitudinal scales A ^trapped). These display a rather strong 
peak in the 2 — 3° range. As noted in the previous section, the transition scale of case 
82126 10:00-12:00 was 1.29°, indicating a rather abrupt transition. However, 37 out of 42 
cases (88% ) have transition latitudinal scales of larger than 2 degrees, which means that 
most cases have more gradual transitions than the examples shown in Olsen et al. [1987]. 




Transition Latitudes 


Figure 6 (bottom panel) displays the transition latitude occurrence histogram. There 
were 31 transitions observed in the southern hemisphere (A tr < 0) and 11 transitions 
observed in the northern hemisphere (A* r > 0). The average transition locations in each 
hemisphere were < A tr >= —5.92 in the southern hemisphere and < A tr >= 8.30 in the 
northern hemisphere. It is also seen that all trapped ion transitions occurred for A > —11° 
and A < 13°, and there is a suggestion of an occurrence peak at about A = —4°. 

2. Correlations the various parameters 

In this section, the correlations among the five parameters and their possible interpreta- 
tion will be presented. These correlations are between: (1) The anisotropy half-width (W) 
and the maximum anisotropy (A max ); (2) The transition latitudinal scale (A,) and the max- 
imum anisotropy (A max ); (3) The transition latitude (A tr ) and the maximum anisotropy 
(A max ); and (4) The penetration ratio (P) for the incoming field-aligned ions and the 
maximum anisotropy (i4 mox ). 


Correlations between Ion Anisotropy Half Width and Maximum Anisotropy 

We might expect that generally the latitudinal half- width W of the trapped ion distri- 
bution should decrease as the maximum anisotropy increases: the more perpendicularly- 

19 



peaked the equatorial pitch angle distribution is, the more confined these ions are in 
mirroring close to the magnetic equator. Figure 7 displays W versus A max for the passes 
considered in the survey. 


In Figure 7, the data show two prominent characteristic relationships: 

(1) cases with very high A m<JI ’s (larger than 7) have lower W’s (less -than 4°). 

( 2 ) cases with very high W’s (larger than 6 °) have lower A maI ’s (less than 5). 

In other words, though there is considerable scatter in the data, the statistical envelope 
of the points is consistent with an inverse relation between W and A m ax- 

The correlation between half maximum anisotropy, W (latitudinal distance between 
locations with A m ax and | A max ), and A max can also be explained as follows. In the model 
discussed above (equation 15), we may write the ratio of the temperature anisotropy at 
location A to its equatorial value as ^ S ^ 

7], ( A )/ 21. ( ? ) V 1 B(\) TIi 9 + /?(A) / & ^ ) 


If we simply let the value oiAx/Amax = 1/2 be associated with a value of ^ L (A)/^ L (eq) = 


7 < 1, we can rearrange this to express the value /2b, 7 as 


Rb, 7 = 1 ~ 


i-1 

2 

Rt.„ — i 


(23) 


where / 2 b, - y is the ratio of the magnetic field strength at the equator to that at the location 
where A\/A max — 1/2. From this last equation we can see that, since 7 < 1 , as tl 
equatorial temperature anisotropy increases, /Zb,- y increases, meaning that the location in 
A is closer to the equator. Thus, this is consistent with our observation of smaller W’s 
associated with higher A’s in Figure 7. In the figure, we have tried three/nttin^curves: a. 
7i/T|| =t 2 \A4maz - 1 (approximation for T± < T d ), b. T±/T\\ = 4v^4 - 3, (approximation 




V © * 


) 


ti ^ 


for Tx. « Td « 60eF), and c. Tj_/T]| = 3vC4. The fitting curves are obtained by finding » y 
the solution for the latitude where >1 = Ai max in equation (23), for a given latitud e^ of / 
ima*i and hence the latitudinal distance between the locations of Ai max a.nd/A max (W) 

As can be seen, curve b. (Tj_ ~ Tj) fits better than curve a. (Tj. T^), rniich may 
indicate that perpendicular ion temperature is close to T^. However, as described before, 
the relation between A and Tj_/Tj| is obtained by crude assumptions: energy dependence 
of the RIMS acceptance solid angle (equation 5) and E m in — 0 (assuming the spacecraft 
is negatively charged or unbiased to significantly simplify the approximation), which will 
affect the correlation between W and A max . For parametric interest, we have tried curve \ 
c (7j_/Tj| = 3 \/A) and found it fits best to the data points. 





Correlations between Ion Latitudinal Transition Scale and Maximum Anisotropy 




V 


cv 

Figure 8 shows that cases with very high A maa .’s have smaller A,’s and cases with 
very large A,’s have lower A max ’s ^while a larger A, indicates a more gradual transitions?^ (Jy 
As in the relationship between A max and W in Figure 7, this is because when the ion 
population has lower A max , the equatorial pitch angle distribution is broad. Hence the 
ion mirror latitudes are spread over a greater latitudinal range which results in a more 
gradual transition. However, the data scatter is rather large in this Figure and the inverse 
correlation is not strong. 


\> 


Correlations between Ion Transition Latitude and Maximum Anisotropy 

Figure 9 shows the tendency that the transition latitude on either side of the equator 
increases with equatorial anisotropy. Note that here the transition latitude is referenced 
from the latitude of A max which is obtained from DE orbit data. Liouville mapping of an 
equatorial bi-Maxwellian distribution peaked in the direction perpendicular to the mag- 
netic field line shows that the distribution becomes increasingly isotropic moving away from 
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the equator. This approach toward isotropy is slower with latitude for lower equatorial 
anisotropy. The correlation between W and A max as shown above shows that the decrease 
of anisotropy with increasing latitude is faster for higher A max cases. On the other hand, 
the larger the anisotropy at the equator, the further it has to decrease overall to attain 

of 1.2. This latter effect is dominant, and the latitude of transition (to 
is larger for larger anisotropjf ||^, as is clearly seen in Figure 9. 

We can also interpret this feature by Liouville mapping. According to our definition, 
the transition latitude is the location where 

j^A=1.2^' Assume the perpendicular temperature is much smaller than 
the detector temperature, we get T±\/T\\\ = 2 \[A — 1 and hence T±\/T\\\ ~ 1.2. Equation 
(22) could be written in the following form: 

A" ^ 

11 ^ & 


Rb> ’ 6 + 6[R T „ - I] 


(24) 


<*> (Ts 


i * 

- -.1 
1 V * 


? X* 
^ «£ «■ 


where is the ratio of the magnetic field strength at the equator to that at the tran- 
sition point. It can be seen that i?£,tr decreases as Rr eq increases, i.e. A tr increases as 
RT' q increases. That is, as the equatorial temperature anisotropy increases, the transition 
latitude will be farther away from the equator. 

Again, according to equation (22), we can fi jyrtltting curve^ for the correlation between 
A tr and A max by assuming different relations between Tj_/Ij| and A max . We have tried two 
different cases, T±/T\\ = 3.2 y/~A and T±/T\\ = 3.3 y/A (the curves vary with the coefficient 
between 2j./Tj| and y/A sensitively.), and found both two curves show the increases of A tr 
as A max increases. The curve with Tx/T\\ = 3.2 y/A fits better in the northern hemisphere, 
and that with T±jT\\ = 3.3 y/A fits better in the southern hemisphere. 

Correlations between Ion Penetration Ratio and Maximum Anisotropy 
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In Figure 10, we plot the penetration ratios iversus equatorial anisotropy parameter 
{A max ). The top panel shows the original points, while the bottom panel shows averages 
within bins of 1 unit of A max ■ As can be seen, there is a clear trend for the penetra- 
tion ratios to decrease with higher equatorial trapped ion anisotropy. This can be readily 
understood by noting that higher equatorial ion temperature janisotropies will generally 
^ 1 / 0> '^ead to greater positive potential values, tending to reflect a larger fraction of the incom- 
ing field-aligned stream back toward the ionosphere, and reduce its penetration into the 
equatorially-trapped ion region. 

It is difficult to find a straightforward correlation between A max and P by utilizing 
Liouville mapping because a variety of factors are involved. For simplicity, we assume 
all the transition points are located at A = 6° such that Rb= 1.06 and further assume 

-ty* 

the parallel ion temperature of cold ionospheric flow to be 0.5eV. For parametric interest, 
in Figure 11, we fix electron temperature at T e = 10eV such that kT e T\\i/e[T e + J||j] varies 
from 0.1 to 5, corresponding to T{= 0.1 to 10 eV. It can be seen, from all the curves, that 
penetration decreases as A max increases and all the data points distribute in the r ange of 
the curves asso ciated with T;=0.1 eV and 21 = 10 eV. -f or ?>•>/ ' 

From the correlations between A max versus A$ and A$ versus P, the trend between 
A max and P may be ^^interpreted as depicted schematically in Figure 12. If we assume 
for the sake of argument that (1) the equatorial electron distributions are approximately 
isotropic or at least uncorrelated with the ion anisotropy, (2) that the actual equatorial 
ion perpendicular temperature is not negatively correlated with the anisotropy, and (3) 
that the effects of the higher energy ring current particles will not change the trend, then 
the higher the anisotropy, the greater the equatorial positive potential peak. The larger 
the positive potential peak, the greater the repulsion and flux reduction of incoming field- 
aligned ion streams. This keeps these field-aligned ions from reaching the equator, and 


results in the decrease of the penetration ratio. 
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4. DISCUSSION AND CONCLUSION 

This report concerns an investigation of core plasma evolution on outer plasmasphere 
flux tubes by data analysis techniques, which we believe has revealed several interesting 
characteristics of the plasma behavior in this inner magnetosphere region. 

In the analysis of DE-l/RIMS core ion data, we have characterized the interface between 
trapped and field-aligned ion distribution regions observed during the unique ’’field-line 
skimming” orbits of the Dynamics Explorer-1 spacecraft. In this analysis, we have charac- 
terized such features as the equatorial ion flux anisotrop y, the penetration of field-aligned 
ionospheric streams into the equatorial region, the scale of the transition into trapped ion 
populations, and the t ransition latitude. We find that there is a trend for the penetration 
ratio, the anisotropy half width, and the transition scale length to decrease with a higher 
equatorial ion anisotropy. 

Liouville mapping of equatorially trapped ions and the reflection of the incoming iono- 
spheric ion streams from the equatorial potential peaks associated with such trapped ions 
have been utilized to interpreted these features. Olsen et al. [1994] used the same mapping 
equations to show observation data points for fitted T± at various latitudes about the equa- 
tor along an L=4.6 field line for pass on January 5, 1984. This distribution of T± agrees 
very well with equation (15) if Tx eq /T\\ eq = 40, and the density agrees with estimated values 
when T± eq /T\\ eq = 18. The agreement of the observed density and temperature distribu- 
tion with the calculated values suggests the distribution of such differential- anisotropy 
produced electrostatic potentials, is likely to be as calculated, with a peak near the equa- 
tor and the decrease toward the ionosphere of 1-2 volts. One of the most important effects 
such potentials will have on the core plasma evolution is the reflection of incoming field- 
aligned ionospheric streams keeping them from penetrating the equatorial region. This 
has been reported by Olsen et al [1987]. The present report provides a statistical survey 
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relating several of the aspects dealt with in the case studies of Olsen et al. [1987,1993]. 

One of the major contributions in this report has 

reveal the effects of the equatorial ion temperature anisotropy on the latitudi- 
nal profiles of plasma parameters along the field line. Similar to simulation results of Lin et 
al. [1994], this study shows results indicating the existence of a positive potential barrier 
around the equator and reflection of field-aligned ion streams. From the interpretation 
of our observational results, we find that as the equatorial ion temperature anisotropy in- 
creases, the equatorial potential increases and hence penetration of field-aligned ion streams 
decreases. 

We may compare the potential drops obtained from simulations with those obtained for 
the same equatorial parameters in the Liouville mapping from equation (16). To compare 
with the simulation results in Lin et al. [1994] (Figure la.) we first assume that the 
equatorial parameters are T e ft: 0.6, Ty\ = 1.1, and = 16. In the simulations, we 
note that the potential drops from the equator to A = 10° and A = 20° are 0.5 and 1 V, 
respectively. Substituting the above equatorial parameters intq^Quatior^^we^would obtain 
potential drops =0.43 at A = 10° and A4? =0.76 at A = 20°. These two sets of values 
are remarkably close. One reason that the potential drops in the simulation are larger is 
that there is an electron temperature gradient which makes an additional contribution to 


the potential^r*^}- 


To make an estimate Jprthe potential drop for a given equatorial ion temperature 
anisotropy observed in the data, we take case 82-126 as an example. According to Olsen 
et al. [1987] the perpendicular ion temperature and A max for case 82-126 are about 8 
eV, so that the equatorial ion temperature anisotropy Rx rq can be estimated to be ft: 5 
from equation 12 (for the case Tj_ eV and A max =10). Hence, the parallel ion 

temperature is about 1.6 eV. If we consideAcases of T„= 0.1, 1, and 10 eV, the estimated 


potential drops at the transition latituda/(A=6.7 0 ) would be 0.02, 0.13 and 0.30 Volts. 
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For these different estimates, with penetration ratio P=0.3, the thermal energies for the 
incoming streams would be 0.017, 0.11 or 0.25 eV (according to equation (20)), respectively. 
The latter estimate is reasonably close to the thermal energies for the ionospheric polar 
wind source. It may also be noted that the streaming ions themselves may alter the 
potential drop from that obtained by simple mapping of the equatorial plasmas alone, 
possibly sharpening the drop and allowing a greater fraction of the total drop between the 
equator and the ionosphere to occur across this transition latitude. 

Here, we summarize the statistical trends in five parameters for 43 latitudinal transition 
events. Indicated trends at this stage include the following: 

1. Transitioi&]£titajityi occur at latitudes below 13°, and about 30% of these occur very 
close to the equator at A < 5°. 

2. A broad range of equatorial ion anisotropies are seen, but the typical values would 
be consistent with bi-Maxwellian temperature ratios of T±/T\ \ of about 3. 

3. The latitudinal scales for the edges of the trapped ion populations display a rather 
strong peak in the 2 — 3° range. 

5. Most cases have penetration ratio P < 1, which implies the ion reflection due to 
positive electrostatic potential layers. P > 1 might indicate an enhancement of the near- 
equatorial flux possibly due to trapping of field-aligned ions in a potential well around the 
equator. These trapped field-aligned ions may not directly be related to the field-aligned 
ions outside the transition latitude. 

6. A higher equatorial A max causes a 
transition scale length and higher transition latitude. 

7. There is a clear trend for the penetration ratios to decrease with higher equatorial 

trapped ion anisotropy. This may be interpreted as follows: larger equatorial temperature 

anisotropies are associated with larger positive electric potential peaks, leading to greater 

repulsion and flux reduction of incoming field-aligned ion streams, hence a decrease of the 
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anisotropy half width and a smaller 



penetration ratio. 
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Figure 1: Spin-time spectrogram of low-energy H + fluxes for a DE-1 pass on May, 6, 1982 
(82-126: Year 82-, Day 126). This case illustrates a rather abrupt transition between 
trapped and field-aligned ion regions and a low penetration ratio of field-aligned flows into 
the equatorial region. 


Figure 2: Spin-time spectrogram of low-energy H + fluxes for a DE-1 pass on March 9, 1984 
(Year 84- Day 068). This case illustrates a more gradual transition and higher penetration 
ratio than case 82-126. 


Figure 3: DE-1 orbit for day 126 of 1982. This orbit shows a rather small L-shell variation 
between A = -10° to A = 0° (A L = 0.13) 


Figure 4: DE-1 orbit for day 295 of 1983. This orbit shows a rather large L-shell variation 
between A = 10° to A = 0° (A L = 0.36). 


Figure 5: Countrates, their ratios and illustration of parameter definitions for the pass 
on day 126 of 1982. (a) shows the time of the pass on maximum flux anisotropy ( A max ) 
and nominal equator(from DE-OA data), (b) shows the penetration ratio of field- aligned 
flows (P), and (c) shows the transition time and transition scale length which indicates the 
sharpness of the transition. 
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Figure 6: Occurrence trends of the ion parameters. 


Figure 7: Correlations between Ion Anisotropy Half Width (W) and the Maximum 
Anisotropy (A max ). The three fitting curves are obtained from equation (23), assuming 
three relations between ion flux anisotropy (A) and ion temperature anisotropy. 


Figure 8: Correlations between ion transition scale (A,) and maximum anisotropy (A moa .). 
Generally, A, tends to decrease as A maz increases, indicating a more abrupt transition for 
a higher A max . 


Figure 9: Correlations between Ion Transition Latitude (A tr ) and Maximum Anisotropy 
( A ma x ) . The two fitting curves are obtained from equation (22). 


Figure 10: Correlations between Ion Penetration Ratio (P) and Maximum Anisotropy 
A max - It can be seen that P tends to decrease as A max increase. 


Figure 11: Fitting curves between P and A ma x f° r different parallel temperatures of the 
equatorial ions. Theses curves are obtained by setting transition latitude At, = 6° and T e = 
10 eV. 
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Figure 12: Schematic representation of increased equatorial trapped ion anisotropy result 
ing in high potential barriers and decreased field- aligned ion stream penetration. 
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Figure 1: Spin-time spectrogram of low-energy H + fluxes for a DE-1 pass on May, 6, 1982 
(82-126: Year 82-, Day 126). This case illustrates a rather abrupt transition between 
trapped and field-aligned ion regions and a low penetration ratio of field-aligned flows into 


the equatorial region. 
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Figure 2: Spin-time spectrogram of low-energy H + fluxes for a DE-1 pass on March 9, 1984 
(Year 84- Day 068). This case illustrates a more gradual transition and higher penetration 
ratio than case 82-126. 
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Figure 3: DE-1 orbit for day 126 of 1982. This orbit shows a rather small L-shell variation 
between A = —10° to A = 0° (AL = 0.13) 
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Figure 4: DE-1 orbit for day 295 of 1983. This orbit shows a rather large L-shell variation 
between A = 10° to A = 0° (A£ = 0.36). 
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Figure 5: Countrates, their ratios and illustration of parameter definitions for the pass 
on day 126 of 1982. (a) shows the time of the pass on maximum flux anisotropy (A max ) 
and nominal equator(from DE-OA data), (b) shows the penetration ratio of field- aligned 
flows (P), and (c) shows the transition time and transition scale length which indicates the 
sharpness of the transition. 
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Figure 6: Occurrence trends of the ion parameters. 
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Figure 7: Correlations between Ion Anisotropy Half Width (W) and the Maximum 

Anisotropy (.4 moa ). The three fitting curves are obtained from equation (23), assuming 
three relations between ion flux anisotropy (A) and ion temperature anisotropy. 
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Figure 8: Correlations between ion transition scale (A,) and maximum anisotropy (A-„ T ). 
Generally, A, tends to decrease as Amax increases, indicating a more abrupt transition for 
a higher A maa . 
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Figure 9: Correlations between Ion Transition Latitude (A tr ) and Maximum Anisotropy 
(A max ). The two fitting curves are obtained from equation (22). 
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Figure 10: Correlations between Ion Penetration Ratio (P) and Maximum Anisotropy 
A max . It can be seen that P tends to decrease as .4 max increase. 








Figure 11: Fitting curves between P and i4 mo * for different parallel temperatures of the 
equatorial ions. Theses curves Me obtained by setting transition latitude A tr = 6° and T e = 


10 eV. 
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Figure 12: Schematic representation of increased equatorial trapped ion anisotropy result- 
ing in high potential barriers and decreased field- aligned ion stre am penetration. 







